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1. Introduction: Feynman's arguments against the variational principle 

The variational principle asserts that for any state in the Hilbert space H of a system with 
Hamiltonian H one finds an energy expectation value that exceeds the ground state energy, i.e. 

£o < ^ ' ' ' 



with Eq the ground-state energy (lowest eigenvalue) of H. If we have a class of variational ansatz 
states |*I'(z)) which are parameterized by a set of parameters z, we can try to find a good approxi- 
mation of the ground state of H within this variational class by finding the parameters z* that mini- 
mize the energy expectation value. The variational method is the basis for a tremendous number of 
highly successful calculational tools in many-body physics. Examples include density functional 
theory [1], Wilson's numerical renormalisation group [2] and the density matrix renormalisation 
group (DMRG) [3]. When applicable, the variational method offers some powerful advantages 
over alternative approaches. It is free of any sign problem that hinders the application of Monte- 
Carlo sampling to many interesting problems, and it is perfectly able to reproduce non-perturbative 
effects and thus easily outperforms perturbation theory. 

In quantum field theory (QFT), however, the variational principle has not met with the same 
success as in other areas of many-body physics. The core reasons for this were identified by Feyn- 
man in one of his last lectures [4]. Feynman pointed out three conceptual issues standing in the way 
of a successful application of the variational principle in (relativistic) QFTs. Feynman phrased his 
first argument as the "sensitivity to high frequencies". This sensitivity is intrinsic to the variational 
method and its attempt to find the "lowest" ground-state energy. While this problem occurs in any 
system containing a large range of interacting energy scales, it is truly catastrophic in relativistic 
field theories. To lowest order, the ground state of a QFT contains the zero-point fluctuations from 
all energy scales and the corresponding ground-state energy is thus dominated by the contribution 
of the high frequencies. In relativistic theories, these UV frequency modes are infinitely abundant 
and create a divergence in the ground state energy (density), which already signals a difficulty that 
the variational principle will face. Any variational approach exploits all variational parameters z 
to obtain the best possible description of the UV degrees of freedom and doesn't care about the 
relatively tiny energy penalty resulting from having an ill-described low frequency behavior. Since 
quantities of physical interest are related to the low-frequency modes, they will typically be very 
badly approximated when using the variationally optimized state Whenever the varia- 

tional parameters affect both the low and high frequencies — unavoidable in interacting theories 
but also true when a real space approach is applied to free theories — this argument can lead to the 
paradoxical situation where the addition of variational parameters provides a worse approximation 
to physical quantities. 

Feynman's second and third arguments concern the lack of suitable variational ansatz states 
allowing high-accuracy computations of observable quantities, rather than an inherent problem of 
the variational method. As a second argument, Feynman remarks that a suitable ansatz for an ex- 
tended quantum field theory should be extensive, i.e. the energy expection value of the Hamiltonian 
with respect to a trial state should be proportional to the volume V of the system. Preferably, 
we would like to work in the thermodynamic limit V — )• oo. For compact systems, one can easily 
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construct a set of variational ansatz states by taking the span of the ground state and a few of the 
excited states of a nearby free Hamiltonian that can be diagonaUzed exactly. For extended systems, 
this approach will fail because the excited states will not be extensive and can thus not contribute 
to the energy expectation value. We thus end up trying to apply the variational method by using 
the ground state of a free theory — a Gaussian state — as variational ansatz, which is equivalent to 
mean field theory. If we do try to devise wave functional which are extensive and non-Gaussian, 
we meet Feynman's third objection: we still have to compute functional integrals in order to cal- 
culate the energy expectation value. This is as difficult as calculating the path integral describing 
the QFT, but in one dimension less, since time does not appear in the Hamiltonian framework. 
For non-Gaussian states, Feynman believed that it is impossible to accurately calculate expectation 
values, as he considered perturbation theory the only means possible to compute these integrals. 
The resulting errors have a strong influence on the optimal state |*F(z*)) obtained by applying the 
variational method and thus on observable quantities derived from it. 

While some progress into non-Gaussian variational methods for relativistic QFTs and gauge 
theories have been developed since Feynman's original lecture [5], we hope that we can contribute 
with a new approach. 

2. A new variational ansatz: continuous matrix product states 

Feynman's second and third arguments ai^e not restricted to relativistic theories, and should 
thus be equally valid for all extended QFTs and even quantum lattice systems. His statement fol- 
lows from his sole consideration of perturbation theory as the only possible method to calculate 
expectation values of non-Gaussian states and expresses a "lack of imagination". Since 1988 [6] 
we have witnessed a stream of results and variational ansatzes for low-dimensional quantum lat- 
tice systems. The most striking development was the introduction of the DMRG [3], which is — 
in retrospect — a variational method within the class of matrix product states (MPS) [7, 8]. The 
basis of the success of the DMRG is the fact that MPS correctly capture the amount of quantum 
entanglement in gapped one-dimensional lattice systems. A better understanding of the behavior of 
entanglement in ground states of short-ranged Hamiltonians — namely the fact that the entangle- 
ment entropy of a region scales as the boundary of that region, with at most logarithmic violations 
for critical systems — has led to the conclusion that ground states of such Hamiltonians live in a 
very small corner of the Hilbert space. This insight was crucial in the development of new varia- 
tional wave functions for strongly-interacting quantum systems that live precisely in this comer and 
have the required scaling of entanglement entropy. The natural generalization of MPS to systems 
satisfying an area law in higher dimensions are the projected entangled pair states (PEPS), whereas 
the logarithmic violations of the boundary law in (1 + 1) dimensional critical systems can be cap- 
tured by the multi-scale entanglement renormalization ansatz (MERA) [9]. These wave functions 
go well beyond Gaussian trial states and allow the accurate and efficient calculation of observable 
quantities. Thus Feynman's objections can akeady be regarded as having been addressed, in the 
case where MPS are applied to relativistic QFTs in conjunction with a lattice regulator [10]. 

These developments have culminated with the introduction of a new variational wave func- 
tional for (1 + 1) -dimensional QFTs: the continuous matrix product state (cMPS) [11, 14]. They 
are obtained as the continuum limit of a certain subclass of MPS and provide a variational class of 
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non-Gaussian wave functionals directly for quantum fields, removing the need for a lattice regu- 
lator. Compelling evidence that cMPS provide a powerful description of the quantum fluctuations 
of quantum fields has been presented in [11, 12]. So far, cMPS have been restricted to the ap- 
proximation of ground states of non-relativistic theories. In this presentation we argue that we can 
apply the cMPS wave functional to relativistic QFTs and illustrate how to overcome Feynman's 
first argument, which is inextricably connected with any attempt to apply the variational approach 
to relativistic QFTs. As for all variational approaches, we work in the Hamiltonian framework and 
explicitly specify the ground state wave functional. While this is not a common approach to tackle 
relativistic QFTs, the Schrodinger formalism for relativistic QFTs is well-established (see [13] and 
references therein). 

A final remark is in order before introducing the cMPS ansatz. In his closing remarks, Feyn- 
man speculated how best to overcome his second reservation and predicted that it should be possible 
to describe a global field state using a reduced set of local parameters. Feynman foresaw the role of 
the density matrix in such a description. It turns out that the density matrix in DMRG has precisely 
the properties envisaged by Feynman: it yields a local parameterization of the global properties of 
a state which is living on the boundaiy of the region of interest. CMPS inherit this as they possess 
a key holographic property: they are pai^ameterized by the (non-equilibrium) dynamics of an aux- 
iliary system — which we could call a boundary field theory — of one lower geometric dimension 
[11, 14]. For (1 + 1) dimensional QFTs, the auxiliary system is zero-dimensional (the Hilbert space 
of the auxiliaiy system can even be chosen to be finite dimensional) and is thus exactly solvable. 
While this is no longer true in higher-dimensional generalizations, we foresee that approximations 
similar to those used in the contraction of PEPSs will still allow for an accurate calculation of 
expectation values. 

While Feynman's arguments ai^e valid both for bosonic and fermionic theories, we focus on 
fermionic theories as these are naturally formulated in terms of creation and annihilation operators, 
which is compatible with the cMPS formalism. We can define the fermionic cMPS class as: 

where [x) are field operators creating fermions of type a at position x with anticommutation 
relations {(//^(x), 1/^3(3;)} = and {Wa{x)^Wp{y)} = 5a,j3 5(x — 3;), Q and /?« ai-e D x D matrices 
acting on the auxiliary system, Traux denotes a partial trace over the auxiliary system, and =^e 
denotes the path ordered exponential. The matrices Q and Ra contain the vaiiational parameters 
and can be position dependent, but we focus on a translational-invariant setting where they are not. 
A derivation of the required algorithmic rules for calculating expectation values of cMPS can be 
found in [11, 14] and we only highlight differences resulting from the anticommutation relations 
of the fermionic field operators. In the relativistic scenario, the two field operators (a = 1,2) 
represent the two components of the Dirac spinor. The state approximates the ground state of 
a relativistic QFT by acting with the field creation operators on the state for which all levels 
are empty (i/^o; |^^) = 0). For free Dirac fermions, the path-ordered exponential should thus fill the 
Dirac sea. 

3. A natural cutoff 

Let's now describe the physical properties of the fermionic cMPS variational class. It is a 
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non-Gaussian class that is both extensive — note the action of the creation operator inside the 
exponential — and allows the exact evaluation of the expectation values of local operators; e.g. we 
obtain (we henceforth use the summation convention on repeated indices): 

for the kinetic energy density, where, in order to obtain real coefficients, we have chosen the con- 
vention = a'' and j8 = a"' for the Dirac matrices. The component vectors (/| and \r) are, 
respectively, the left- and right-eigenvectors of the transfer matrix T = 1+ l0Q + Ra®Ra, 
corresponding to eigenvalue zero [11, 14]. We focus on the kinetic energy density as it is the domi- 
nant term in the UV region, which is the region responsible for divergences and for Feynman's first 
criticism. As long as the D x D matrices Q and Ra have finite entries this expression will be finite 
and is thus regularized. 

A better understanding of this regularization is gained by looking at the momentum occupa- 
tion in a cMPS: {x\Vai!<^)Vpi^')\x) = ^(^~^')"a,/3(^) [1^]. The momentum occupation number 
"a,j8 is the Fourier transform of Cap (x), where 

Ca.p{x) = d{-x){l\{l^ Ra^^ {Rp(^l)\r)+e{x){l\{Rp^iy^{l(g>Ra)\r) 

and d{x) the Heaviside function and T = Q(i$ 1+ Q — Ra ® Ra, where the last minus sign 
originates from the Fermi statistics of the particles. There will not be any disconnected contribution, 
as we require {x\¥a\z) = 0- The behavior of na,j} (k) for large k is determined by the continuity and 
differentiability of Ca.p, in particular around x = 0, which is the only point where differentiability 
of the expression above is not trivially guaranteed. Since Ca.pix) is a continuous function, its 
Fourier transform decays as nap{k) < t&{k^^) for |^| — oo. Continuity of the derivative of Ca.p{x) 
at a; = requires {l\ {Rp,Ry} ® [RajRy} k) =0 (Va,j8), which is satisfied by choosing all matrices 
Ra nilpotent and anticommuting. The second derivative of Ca.p (x) at a; = is then automatically 
continuous, from which one can conclude that na^p{k) < ^(^^'*) for |^| — oo. While a faster 
decrease of the momentum occupation number imposes additional constraints on the matrices Q 
and Ra the cuiTcnt behavior already ensures a finite kinetic energy. The region in momentum space 
where the fc^^ decay behavior sets in defines a soft momentum cutoff A. 

4. Curing the sensitivity to high frequencies 

We can now investigate how Feynman's first objection manifests itself for the cMPS ansatz. 
The problem is situated in a cMPS's ability to describe a scale transformation x cx (c > 0) 
by an equivalent transformation Q' = cQ and = y/cRa- Since this transformation does not 
change (/| and |r), the kinetic energy per unit length will simply be multiplied by a factor c^. In 
renormalizable theories, the kinetic energy has the highest scaling dimension, together with other 
terms with dimensionless coupling constants. These are thus the dominant terms in the UV region. 
However, in contrast to the non-relativistic case, the relativistic kinetic energy is not a positive 
definite operator and can acquire a negative expectation value. \i\x) is a cMPS for which the sum 
of terms with highest scaling dimension has a negative energy expectation value, then the total 
ground state energy can always be decreased by a scale transformation with c sufficiently large. 
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low frequencies 




momentum k 

Figure 1: Hypothetical momentum distribution of an optimal cMPS for a free fermionic theory: high- 
frequency degrees of freedom are well-approximated up to a cutoff A, after which the momentum occupation 
decays as k^"^. Also shown is the effect of a scale transformation. 

Our variational method will thus try to push c — )• oo, in order to approximate the divergent (kinetic) 
energy of the exact solution. Under such a transformation, the momentum occupation changes to 
^'a jS ~ i^/^) the intrinsic cutoff determined by n' is given by A' = cA. 

This change of scale will be accompanied by a worse description of the low frequency region, 
as predicted by Feynman. The precise underlying cause for this effect in our variational class is 
that a cMPS can only accurately describe states with a finite amount of entanglement. The maxi- 
mal entanglement entropy in a one-dimensional system with energy gap A and energy cutoff A will 
roughly be given by S ~ log (A/ A), and a cMPS with D proportional to ^(exp(5)) should suffice 
to provide a good description [16]. If D is too low to obtain a good approximation of the exact 
ground state, the variational method wiU make compromises in that part of the frequency spectrum 
that contributes least to the ground state energy, i.e. the low-frequency region. In non-relativistic 
systems, the cutoff is set by the particle density or thus by the chemical potential. But in a relativis- 
tic Hamiltonian, there is no physical cutoff and we only have the intrinsic momentum cutoff A of 
the cMPS. If we start from a cMPS with negative energy expectation value, the variational method 
can quickly lower the energy by shifting the cutoff to A' = cA with c — oo. As c goes to infinity, all 
low-energy modes will eventually fall into the region that is poorly described and the description 
of any observable quantity will be completely wrong for every finite value of the bond dimension 
D. This is schematically illustrated in Fig. 1. 

A solution is now straightforward as we can prevent c from running to infinity by imposing 
a constraint on the matrices Q and Ra'- since Q has the dimension of momentum and Ra has 
the dimension of the square root of a momentum, constraining the nonn of Q and Ra prevents 
c from running and regularizes the resulting theory by introducing a scale, i.e. a dimensionful 
parameter, into the system, similar to what happens in analytical regularization techniques or lattice 
regularization. In the sequel, we will constrain the norm of the commutator [2,/?o;] by fixing the 
expectation value of (di//^/djc)(di///dx) [17]. Hereto, we add this term to the Hamiltonian with 
a Lagrange multiplier 1/A, i.e. ^cutoff = Ax{A'^'^ /dx){d'^/Ax). This apparently arbitrary 
choice is motivated by the requirement that the constraint needs to penalize high values of the 
momentum k, to which [2,/?o;] is related by the calculational rules of cMPS. //cutoff will give a 
contribution in momentum space, which is low enough to ensure a finite result in combination 
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with a momentum occupation that decays as k^'^. It is, however, strong enough to penahze high 
frequency modes, even the ones that give a contribution — to the (kinetic) energy. Put differently, 
it is a positive definite term with a higher scaling dimension than the relativistic kinetic energy. 
As such, it is non-renormalizable, which by means of the renormalization group indicates that it 
will be irrelevant for the description of the low-frequency modes and cannot strongly influence 
the expectation value of observable quantities. Note that it does respect the chiral symmetry of 
the kinetic energy term. It does of course break relativistic invariance, which is inevitable when 
introducing a momentum cutoff in a Hamiltonian framework. We expect that any other norm 
constraint with similar properties and respecting the symmetries of the system should also work. 



5. Application 1: free Dirac fermions 



We now illustrate our arguments by ap- 
plying them to relativistic fermion models. 
As a benchmark, we first consider free Dirac 
fermions with Hamiltonian density 



f i .t/ N vdV>^, , i diif' 

'^-=-2'^'W^^d^W + 2^ 
+ m^r\x)o^^r{x), 




with Dirac matrices chosen as described 
above, and m the fermion mass. In the exact 
ground state of Hd + /^cutoff > this term will 
actually introduce a sharp cutoff at ^cutoff = 
A(l/2 + (1/4 + mVA2)i/2)i/2^ which is 

equal to A up to corrections of ^(ni^/A^). 
The cMPS ansatz will not be able to repro- 
duce this sharp cutoff because it decays as 
k^'^. Indeed, this cutoff is not expected to 
be reproduced very well, because the new 
hamiltonian is gapless at /: = ib^cutoff- How- 
ever, this is not a problem, as we do not 
expect these high-frequency modes to influ- 
ence physical properties. Note that both ze- 
ros in the dispersion relation occur at physi- 
cally different momenta and do thus not result in fermion doubling. 

Since we do not aim at reproducing the exact solution in the high-frequency regime, we can 
not compare the corresponding energy as a measure of the accuracy of our solution. Instead, we 
have calculated the momentum occupation of the exact positive (particle) and negative (antiparticle) 
levels according to the definitions {a^{k)a{k')) = 5{k' -k)n++{k), {P {k)b{k')) = 5{k' -k)n— (k), 
{a^ {k)b{k')) = 5{k' — k)n^^ {k), with d {b) the annihilator (creator) of particles (antiparticles). 
The results corresponding to the optimal cMPS are shown in Fig. 2. The exact solution has the 
Dirac sea filled {n [k) = 1) all the way up to ^cutoff> after which n [k) = for |^| > fccutoff> 



Figure 2: Momentum occupation of the antiparticle 
levels n (k), the particle levels n++(A;) and the mix- 
ing in a cMPS approximation of the Dirac 
field with mass m. The auxiliary space of the cMPS is 
C'^ ®<C^ ®£P . The vertical Une indicates the position 
of the exact cutoff. 
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and {k) = ri^^ {k) = 0, V^. These results were obtained using the cMPS ansatz where Q and 
Ra act on an auxiliary Hilbert space (g) C^, where the first two two-dimensional Hilbert 
spaces accommodate auxiliary fermions which are used to impose the anticommutation relations 
on Ra- It is clear from these results that the low-energy behavior is approximated very well for 
the massive Dirac theory, and the accuracy greatly increases by increasing D. As anticipated, the 
cutoff behavior is approximated less well. In the case m = the theory is critical and the low- 
energy behavior is also approximated less well. This result is familiar from MPS solutions for 
gapless lattice models. However, from the fact that |«+^ | for m = 0, we see that the algorithm 
automatically converges to a cMPS respecting chiral symmetry, except at D = 2. 

To give a non-trivial example of what 
can be done with our approach we have also 
calculated the Casimir energy of the massive 
Dirac field. We simply recycle the matrices 
Q and Ra from the simulation above, and 
add suitable operators B{x) to the ansatz at 
the location of the 'plates' or defects (x = 
and X = L), which impose the correct bound- 
ary conditions. The boundary conditions will 
only fix a part of these operators, the remain- 
ing elements can be used to optimize the en- 
ergy with fixed Q and Ra- AH boundary ef- 
fects can be incorporated in the boundary op- 
erators B- We used the typical bag-model 
boundary conditions [18]. In Fig. 3 we show 
the energy density for a particular- configu- 
ration of 'plates' in the (1 + 1) dimensional 
free-fermion model. A clear manifestation 
of the Fermi surface at a finite momentum 
^cutoff is present in the form of Friedel oscil- 
lations. The Casimir energy Eq as a function 
of the distance L between the plates is plot- 
ted in Fig. 4. The presence of the momen- 
tum cutoff, and thus of the finite particle den- 
sity, also introduces a strong oscillatory be- 
havior in Ec{L), which was already observed 
in studies of the interaction energy between 
defects in one-dimensional quantum liquids 
[19]. Local minima correspond to values 
of L where the number of allowed modes is 
such that the density of fermions between the 

plates is exactly equal to the density of fermions outside the plates. The sharp maxima appear when 
this condition is most strongly violated. In the limit /ccutoff °°, the density of fermions is infinite, 
both in between and outside the plates, and the equal density condition is always satisfied. So the 



Figure 3: Energy density e for Dirac fermions with 
m/A = 1/10 in a system where 'plates' are present at 
position X = and x = 50/A. These plates enforce 
the bag-model boundary conditions. The ground state 
energy density eo in the infinite vacuum is plotted for 
comparison. 



0- 



-0.05 



-0.05 




20 



40 



60 



80 



100 



L A 



Figure 4: The total Casimir energy as a function of 
the distance L between the 'plates'. For m = 0, the 
exact Casimir energy, both in our model with cutoff, 
and analytically through zeta function regularization, 
is also displayed. 
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physical Casimir energy, which is expected to be cutoff independent, can be found by the envelope 
of the local minima. This is illustrated for the m = case, where the exact solution in our model 
is compared to the value of the Casimir energy for massless Dirac fermions in ( 1 + 1 ) dimensions, 
as calculated with zeta-function regularization, i.e. Ec{L) = —j^- Note that the Casimir energy 
will always have an asymptotic exponential decay in cMPS calculations, but that it can be well 
approximated at intermediate values of L. 

Because the Casimir energy is a difference of energies approximate results can be lower than 
the exact solutions. This is clearly the case in Fig. 4. We attribute this effect to the fact that the 
additional degrees of freedom present in the boundary vectors allow one to further optimize the 
energy in their immediate vicinity. The qualitative behavior of the energy is already reproduced by 
this simple approach. 

6. Application 2: the Gross-Neveu model 

As a final proof of principle, we study a theory with interactions. One of the most important 
models for one-dimensional relativistic fermions is the Gross-Neveu model, as it shares many fea- 
tures with QCD [20], including, asymptotic freedom and spontaneous breaking of chiral symmetry. 
The hamiltonian density for the A/^-flavor Gross-Neveu model is given by 



where the ;c-dependence of the field operators has been omitted for brevity and there is an implied 
summation over the flavor index a = l,2,...,N. One must not forget to apply normal ordering 
when deriving an interacting Hamiltonian from the relativistic path integral, which is a coherent- 
state path integral for fermionic theories. As a variational ansatz we employ a product state of cMPS 
states across the different fermion flavors. Because the exact ground state has flavor symmetry 
(and actually 0{2N) symmetry), the nearest product state should also be invariant under [21]. 
We can thus use the same cMPS for every flavor. This amounts to a Hartree-Fock approximation of 
the theory, where the self-interaction of the flavor is treated exactly, and the self-consistent mean- 
field approach is only applied to the interactions between different flavors. We add the same cutoff 
term ^cutoff to the Hamiltonian for every fermion flavor, so as to respect the flavor symmetry. Since 
this term introduces our regularization parameter A, we know that the coupling constant g will have 
to depend on A in order to have a consistent theory. In the N limit, we can solve this problem 
exactly, and we obtain the well-known result for a = {xlw'^ (^~'¥\x) 



where /ccutoff ~ A if |Aa| ^ A. This indicates that the cutoff fixing term ^cutoff has no effect other 
then what it is meant to do, i.e. introducing a cutoff. With the current Hartree-Fock ansatz based 
on cMPS, we can calculate an approximation for any A and N. In principle, we can describe 
the complete wave function of all N flavors with a single cMPS, but this is computationally more 
demanding as the dimension of the auxiliary space needs to grow exponentially with the number of 
flavors. Numerical results with the mean-field approach are illustrated in Fig. 5. At strong coupling 
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Figure 5: Expectation value of d = {x\ i/^ff'i// \x) in the Gross-Neveu model as function of A(A) forN - 
A fit of the form cie^'-/^ to the numerical results for A^^ < 1 at D = 16 results in C2 
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-0.072' 
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to be compared to the exact values ci —2 and cj — Ti (see main text). 



(A > 1) they agree very well with the exact result. The discrepancies between the exact solution 
and the cMPS approximation for = oo are clearly finite-D effects. They become more pronounced 
as A a/ A gets smaller, since A a is exactly equal to the mass gap in the N = oo limit. 

7. Conclusion 

We have developed an extension of the cMPS variational class appropriate for fermionic 
(1 + 1) -dimensional relativistic field theories. Since cMPS have a built-in cutoff, they offer a new 
way to regularize quantum field theories. We have explained how to ensure that cMPS do not suffer 
from Feynman's objections concerning the application of the variational principle. Additionally, 
our approach is free from fermion doubling and sign problems. We have demonstrated the appli- 
cability of our variational approach by reproducing the known results for free Dirac fermions and 
provided two nontrivial applications of our method to the Gross-Neveu model, where we observe 
chiral symmetry breaking and absolute scaling of the chiral parameter, and the Casimir effect, 
where we are able to reproduce the qualitative features of the Casimir energy. 
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